Reflectivity Anisotropy Spectra of Cu- and Ag- (110) surfaces from ab initio theory 
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We present the surface Reflectivity Anisotropy spectra of the relaxed (lxl) (110) surface of Cu 
and Ag as obtained by ab initio calculations. We are able to disentagle the effects of the intraband 
and interband parts of the bulk dielectric function on the bare dielectric anisotropy of the surface. 
We show how the position, sign and amplitude of the structures observed in such spectra depend 
on the above quantities. The lineshape of all the calculated structures agree very well with the 
ones observed experimentally for samples treated by suitable surface cleaning. In particular, we 
reproduce the observed single peak structure of Ag at high energy, found to represent a state of 
' ^ ■ the clean surface different from the one giving the originally observed double peak structure. This 

C/3 ' results is not reproduced by the 'local field' model. 



The Reflectivity Anisotropy Spectroscopy (RAS) is a powerful technique to probe the optical properties of an 
anisotropic medium like the surface of an otherwise isotropic bulk crystal. The RAS is measured as the relative 
difference — of the complex reflectivity r along two perpendicular axes on the surface. Despite its vast application 
in semiconductors surface science since the pioneering work by Aspnestl, the first RAS spectrum of a metallic surface 
appeared as late as in 1993 for the Ag (110) surfacea. The characteristic resonance-like structure at the bulk plasma 
frequency, i.e. at uj p = ~ 3.8 eV, has been considered a fingerprint of the Ag-(llO) surface. More recently, the research 
on the optical properties of metallic surfaces has gained renewed impulse and in the las.t.±hree years the (110) surface 
O \ of Ag and Cu have been investigated both by optical and photoemission spectroscopyaQ. 

The RAS results essentially agree in establishing the existence of one peak at low energy (< 2.5 eV) of similar 
y—{ \ shape in both Cu and Ag and of a further structure at higher energy of quite different shape. In particular, some 
i authors point out thatpthe lineshape of the high-energy peaks depends critically upon the surface treatment during 
0^ ' the sample preparationHQ. This experimental aspect is very important for a correct comparison of RAS spectra with 
, theoretical models and calculations. 

A major problem with the interpretation of the RAS spectra is that the optical anisotropy introduced by the surface 
is difficult to extract from the measured reflectivities without a solid quantitative theoretical basis. An attempt to do 
so has been made from a phenomenological point of viewQ. Experimentally one can observe the evolution of a given 
feature in the spectra from the clean to the gas-covered surface. In fact coverage by, e.g. oxygen, is known to wash 
i out progressively optical structures originating from surface-to-surface state transitions, as a function of time and/or 
0$ ' amount of oxygen. 

RAS structures observed in the noble metals are currently interpreted in terms of interband transitions between 
surface or bulk states by speculating on available results of existing band structure obtained by folding bulk bands or 
by interpolation methods. Another tool of interpretation is the so-called 'swiss cheese' modeB, describing the optical 
response of the surface by layer-depending dipolar interactions. This phenomenological model partially succeeds, 
at least qualitatively, to explain the RAS peaks observed in Ag aad Cu under suitable experimental conditionsflfJ. 
The lineshape, however, depends on the parameters of the modefl. Only, very recently first principle calculations 
of the optical properties of noble metals surfaces have started to appeanljO, opening the possibility of a direct 
quantitative interpretation of the measured optical spectra. Such self consistent calculations yield, on the same 
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footing, the electronic structure and the dielectric functions e^fgi) e [ooi] along the directions [110] and [001] of the 
(110) surface . The dielectric function is computed from interband transitions among electronic states. From these, 
Ae= e^jQi —£[ooi]: known as Surface Dielectric Anisotropy (SDA), is derived. So far, comparison of the calculated SDA 
with direct experimental RAS results is made through a reliable semi-phenomenological modelBEl that transforms 
the RAS spectrum into the SDA via the bulk dielectric function e b normally obtained by independent measurements. 
Attempts to interpret the RAS spectrum in terms of bulk or surface electronic transitions from the SDA structures 
may, however, be incorrecuij. In fact, the passage from the anisotropic dielectric function to the RAS spectrum 
involves bulk quantities that re-shape much of the bare metallic SDA spectrum. At low energies, as shown below, 
RAS and SDA have different signs due to intraband transitions, at variance with the case of semiconductors. This 
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point is even more important in view of recent controversial results for the sign of the RASM3. Therefore, as we 
aim to show in this Letter, only the calculation of the RAS spectrum from consistently computed quantities allows 
to distinguish between RAS structures arising from the bare SDA from those related to bulk effects. 



We shortly summarize the fundamental analytical expressions necessary for the analysis of a RAS spectrum. The 
surface reflectivity anisotropy — , written for the present surface geometry, is — = 2 11101 10011 . The RAS spectrum 

r r r [no] +r [ooi] 

, as obtained directly from reflectivity measurements, is usually given as the real part of the above complex quantity, 
namely: 

R f Ar,_ Ae 2 (e\-1) e|Ae x 

^ °^-l) 2 + (4)2 (e 6_ 1)2 + (e 6 )2 i' 

where subscripts 1 , 2 refer to the real and imaginary parts, respectively and C is a constant that is of no importance 
here. 

In metals, all the e's in the formula above contain terms due to interband as well intraband transitions. The latter 
ones give rise to the so-called Drude tail dominating at low energies. We assume that the anisotropy of the intraband 
part of Ae is negligible, which allows us to neglect intraband effects tout court in this term. In fact there is, so far, 
some experimental evidenced that the intraband anisotropy of the surface affects the spectrum of Cu only up to 1.5 
eV. The inclusion of intraband effects in e b is, instead, of paramount importance, as we will show below. Therefore our 
calculations contain the ab initio computation of the interband part of all the surface and bulk e's plus the inclusion 
of the Drude tail in|-e, b computed by the standard relations involving the plasma frequency cu p and the relaxation time 
r of the bulk metaO. 

The computation of ab initio quantities, i.e., electronic properties and intraband part of the dielectric fitactions, has 
been carried out self-consistenly by the Full Potential versionEj of the Linear Muffin Tin Orbital methodlij within the 
Local Density Approximation!!!! for the exchange and correlation potential!!]. The surface geometry is obtained within 
a repeated slab scheme including 11 atomic plus 6 vacuum layers. The reciprocal space integrations are performed 
with the analytical tetrahedron methodlj. Meshes as dense as 752 and 256 k-points in the irreducible bulk and 
surface Brillouin Zone, respectively, give converged values for the dielectric functions calculated as matrix elements of 
the momentum operator among occupied and unoccupied electronic states. The relaxcdEEl, unreconstructed geometry 
of the Cu and Ag (110) surfacestil has been taken into account. A more complete description of the details of the 
calculations will be published in a future work. 

In Fig. [jjand Fig. ||, we have plotted the calculated RAS, (Re{^-}), and SDA, (e^ygj — £[ooi])i f° r Cu and Ag, 

respectively. To disentangle the effect of the intraband from interband part of e b in the RAS we have also plotted in 
the upper panel of Fig. |l| and Fig. || Re{ ^r} calculated without the intraband (Drude) term in e b . In Fig. || we show 
the plots of the calculated real and imaginary part of e b , separated into total and interband parts, for bulk Cu (upper 
panel) and Ag (lower panel). The plots shown in these figures illuminate very clearly that in both metals, the overall 
effect of e b is to re-shape profoundly the bare dielectric anisotropy of the surface Ae, depending on the energy and 
on the peculiarities of the bulk dielectric function of each metal. The RAS spectra of the two metals may be divided 
into a low (< 2.5 eV) and high (between 3 and 5 eV) energy part, according to the prominent spectral features. The 
low energy peak has a similar shape in both metals whereas the high energy structure looks quite different. 

Let us now focus on the two metals separately in order to highlight how characteristic quantities of the bulk, e.g. 
the plasma frequency, affect the RAS spectrum. In Fig. [j] (upper panel) the RAS spectrum of Cu, with and without 
the Drude term in e b may be compared with the SDA in the lower panel. From this one sees very clearly that the 
low energy peak in the SDA undergoes a sign reversal and a little energy shift due to the intraband part of e b . This 
effect decreases with increasing energy until the interband part of e b takes over. For energies larger than ~5 eV the 
effect of the Drude tail dies off almost completely. The double peak shape of the high energy structure is essentially 
determined by the interband part of e fc between 4 and 5 eV (see Fig. ||, upper panel) whereas the Drude part of e\ 
amplifies the negative part of this peak between 4 and 4.5 eV. 

In Ag the effect of e b is even more dramatic, as can be seen from the two curves in the upper panel of Fig. by 
comparing them with the SDA plot in the lower panel. The low energy part of the spectrum (see the inset of Fig. g) 
is determined by the same interplay of bulk intraband e b and SDA as observed in Cu. The RAS spectrum is, however, 
dominated by the effect of the bulk plasmon producing the overwhelming negative structure at about 3.2 eV. This is 
due to the steep behaviour of e b close to the treshold of the onset of the interband transitions, as can be inferred from 
Fig. H (lower panel). Afterwards, the spectrum is only sligtly affected by the Drude term, as in the case of Cu. 
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Let us now come to the direct comparison of our. calculated RAS spectra with the corresponding experimental ones. 
As pointed out very clearly by extensive studiesa'tj, the lineshape of the observed structures may vary drastically not 
only upon epxosure in air but also upon the number of cleaning cycles of the surface. p n 

RAS spectra of Cu and Ag clean (lxl) (110) surface have been measured by various groupsOfllJO in different 
experimental conditions. All these spectra display the low energy peak at ~ 2.1 and 1.8, in Cu and Ag, respectively 
and the high energy structure located between 4 and 5 eV, in agreement with our findings above. However, the shape 
of the latter structure depends upon treatment of the surface. In particular, a recent paper on the effect of sputtering 
on the RAS spectrum of the (110) surface of Cu shows that the shape of the double-peak structure at high energy (see 
Fig. [l] ) varies very much with sputtering, holding as fixed the 'node' of the structure (see Fig.l of Ref.tJ). This node 
occurs also in our calculations at essentially the same energy. We emphasize that the position and overall shape, af 
the two structures in the RAS of Cu presented in this letter agree quite satisfactorily with the experimental resultsHfl. 

Also for the Ag clean (lxl) (110) surface our calculated RAS spectrum in Fig || gives satisfactory agreement with 
the experimental findings apart from the shift of about 0.5 eV toward lower frequency in the location of the high energy 
structure. We attribute this discrepancy to the value of the threshold of the interband transitions that amounts to 
3.3 |©y in our calculations, as can be seen in Fig. || (lower panel). This value is 0.5 eV smaller than the experimental 
ondi-3. tEhis discrepancy has been found also in another calculation for bulk Ag, carried out using the pseudopotential 
methodE3; it must be due to the shortcoming of Density Functional Theory, that sometimes underestimates transition 
energies.) The lineshape of this structure is surprisingly in agreement with the one obtained in very intensively 
polished samples, as discussed abovetJ. The double peak, resonance-like structure measured in freshly polished Ag 
surface, instead, does not occur in our ab initio calculations, since they.-.do not include the long-ranged local-field 
effect. This is described by the 'swiss-cheese' phenomenological modela where the noble metal is modeled by a 
composite medium made by a jellium-like background (the sp electrons) screening the embedded charged spheres, 
centred about ions carrying dipoles (due to the d electrons) interacting among themselves. The effect of the reduced 
symmetry of the surface on the dipole-dipole interaction leads to different plasma frequencies for the different light 
polarizations. This is the origin of the derivative-like lineshape observed in Ag RAS at high energy. Its absence at 
surfaces treated by several sputtering cycles may be understood by the introduction of defects {e.g. vacancy islands) 
by the treatment itselfu. The geometrical roughening of the surface may destroy the long range order, suppressing 
the effect of dipolar interactions and leaving only the dominant effect of the local environment, well described by our 
one-electron calculations^. 

To conclude and summarize, we have in this Letter calculated ab initio the RAS spectra of Cu and Ag (110) surfaces, 
showing which structures are due to the surface dielectric anisotropy and which are due to bulk effects. Moreover we 
have disentagled intraband from interband contributions of the bulk dielectric functionJn the RAS spectrum, showing 
that it is incorrect to interpret RAS_suectra from the knowledge of the SDA onlytil. For the first time the RAS 
structures observed at metal surfacesaQ have been accurately reproduced by ab initio calculations. The measured 
lineshapes in the Cu and Ag RAS depend critically on the conditions of the surface. As discussed exstensively in 
some papers, sputtering, annealing or exposu£|£,in air may change not only the amplitude but also the shape of the 
structures as found ipr the high-energy peaksErQ in Cu and Aa. n We conclude that our results and those obtained by 
the local-field modeB refer to samples whose surface conditionsEm emphasize short-range and long-range arrangements 
of the atoms on the surface, respectively. 
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FIGURE CAPTIONS 



FIG. 1. RAS (upper panel, solid line) and SDA (lower panel) spectra for the Cu-(llO) surface. The dotted curve in the 
upper panel is the RAS obtained excluding the Drude tail from e . 



FIG. 2. The same plots as in Fig. m for the Ag-(llO) surface. The inset shows the low energy peak on a magnified scale 



FIG. 3. Calculated Real and Imaginary parts of e b for Cu (upper panel) and Ag (lower panel): total (bold full line) and 
interband (dotted line) e\, total (dashed line) and interband (dash-dotted line) e\. The bulk free-electron lo p and relaxation 
time r are taken from Ref. 11 to be 9.17 eV, 6.9 10 -15 sec -1 and 9.22 eV, 31 10 -15 sec -1 for Cu and Ag, respectively. 
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